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ABSTRACT 

We estimate the thermal conductivity of a weakly coUisional magnetized plasma with chaotic magnetic 
field fluctuations. When the fluctuation spectrum extends over two or more decades in wave-vector, we 
find that thermal conduction is very efficient; the conduction coefficient is only a factor ~ 5 below the 
classical Spitzer estimate. We suggest that conduction could play a significant role in cooling flows in 
clusters of galaxies. 

Subject headings: galaxies: clusters: general — cooling flows — conduction — magnetic fields 



1. INTRODUCTION 

H ot X-ray-emit ting gas is ubiquitous in clusters of galax- 
ies (Barazin 1988). Since the X-ray emission is energeti- 
cally important in many clusters, it is believed that a sig- 
nificant amount of mass must co ntinuously co ol and drop 
out of the intracluster medium (Fabian 1994). The mass 
deposition rate is estimated to be as m uch as several hun- 
dred M(7)yr~^ in some clusters (e.g., David et al. 2001 



A llen et al. 200l| ). 

Direc t evidence fo r the cooling gas has, however, been 
scarce ( Fabian 19941). In particular, recent observations 
with XMM- Newton ( (Bohringer et al. 2001|; [Molcndi fc Pizj 
zolato 2001) and Chandra (Fabian et al. 2001) have failed 



to find the multi-temperature gas one expects in a cooling 
flow. The observations suggest that mass dropout may be 
less significant than previously thought. A reduced level 
of mass dropout is possible if there is a source of heat to 
replace the energy that is lost through X-ray emission, but 



no cl ear heat source has yet been identified (Fabian et al 



2001). 



The inner region of a cluster [R ~ few x 10 kpc) where 
mass dropout appears to be occurring is typically cooler 
than the rest of the cluster. Therefore, an often-discussed 
source of heat is thermal conducti on from the hot outer 
regions of the cluster to the center (|Binncy fc Cowic 1981 



Tucker fc Rosncr 1983|; |Bertschingcr fc Mciksin 1986|; |Breg 



man fc David 1988|; |Gaetz 1989|; |Rosncr fc Tucker 1989 
Pistinncr fc Shaviv 1996|; [Santos 200 ip. While the idea 



is attractive, it requires extremely efficient conduction, 
which is considered p roblematic. 

In a classic paper, ^pitzer (1962| ) showed that thermal 
conduction in an unmagnetized plasma has a diffusion con- 
stant^ Ksp ^ AVicoui = Aut - 4x10^^ nZl cm^ s"\ 
where tcoui — X/vt is the mean free time between Coulomb 
collisions, and A and Vt ar e the mean free path and the 
thermal speed of electrons (Cowic & McKee 1977; Ettorj 
fc Fabian 2000| ): A ^ 30 T^^^ 
s-\ Here, Ti : 



4 X 10« T, 



9 



kpc, Vt ~ {kT/mef^ ~ 
= A:T/10 keV is the scaled 



temperature, n_3 = n/10^'^ cm^'^ is the scaled electron 
number density, and we have used an average value for 
the Coulomb logarithm. In A ~ 38. 

The time required for heat to diffuse conductively across 
a radius R = 100i?2 kpc is given by isp ~ -R^/^Sp ~ 



10^ ^'^n. 



3R2 yr. For conduction to have a signif- 
icant effect on a cooling fiow, the conduction time must 
be comparable to the cooling time tcooi of the gas. Ta- 
blc 1 lists represen tative data f or two clusters. H ydra A 
dDavid et al. 2001] ) and 3C295 ( [Allen et al. 2001| ), at two 
characteristic radii, 100 kpc and 10 kpc. Columns 6 and 
7 give icooi and tsp. We see that if thermal conduction 
in a cluster is as efficient as in Spitzer's theory, or even if 
it is a factor of a few less efficient, heat conduction will 
have a strong effect on the energetics of a cooling flow and 
perhaps will shut off mass dropout. The main problem 
with this idea is that the gas in a cluster is likely to be 
magnetized, and conventional wisdom says that magnetic 
fields severely suppress conduction relative to the Spitzer 
level. This is the topic of the present Letter. 

We discuss in §2.1 the theory of conduction in a tangled 
magnetic field as developed by Rechester fc Rosenbluth 
(1978, RR) and Chandran fc Cowley (1998, CC); the the- 
ory predicts that the coefficient of thermal conduction is a 
factor ^ 100 — 1000 lower than the Spitzer coefficient. We 
then present in §2.2 and §2.3 an extension of the theory to 
a turbulent medium; we show that if turbulence extends 
over a factor of 100 or more in length scale, thermal con- 
duction is almost as efficient as in Spitzer's theory. We 
conclude with a brief discussion in §3. 



THEORY OF THERMAL CONDUCTION IN A WEAKLY 
COLLISIONAL MAGNETIZED GAS 



2.1. 



Conduction in a Chaotic Magnetic Field With a 
Single Scale 

In the presence of an ordered magnetic field, conduc- 
tion is anisotropic. Electrons stream freely parallel to 
the field line, so the parallel diffusion constant is almost 
equal to the Spitzer value: K|| ~ Ksp/3. (The factor of 
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1/3 is because diffusion is in one dimension rather than 
three, see CC). Perpendicular to the field, however, elec- 
trons follow circular Larmor orbits with radius pe <S^ A. 
Since an electron moves only a distance ~ pe in each scat- 
tering, the perpendicular diffusion constant is given by 
pl/tconi ~ (Pe/A)^KSp < Ksp- For a galaxy cluster 
with a magnetic field of ~ 10^^ G, we have pe ^ lO^^^A, 
and K±_ is effectively zero. 

Thermal conduction behaves very differently when the 
magnetic field is chaotic. The theory for a tangled field 
with a single coherence length Ib was developed by RR 
and has been recently revived in the astrophysical context 
by CC. Since the field is chaotic, the separation r of two 
nearby field lines must have a Lyapunov-like scaling as a 
function of distance I along the field: 



r ~ To exp(//LLyap), 



(1) 



where ro is the initial separation of the two lines. As there 
is only one characteristic scale in the problem, namely Is, 
we expect ^Lyap ~ Ib- 

Following RR and CC, let us consider the evolution of 
a compact cloud of electrons of initial size pe- With time, 
the electrons diffuse parallel to the field, with a diffusion 
constant Ky . As the electron cloud spreads out, its perpen- 
dicular extent diverges exponentially according to equa- 
tion (|l|). Thus, by the time the electrons have diffused 
a Rechester-Rosenbluth distance L_rh ~ Ib^^{Ib / Pe) ^ 
iOls along the field line, their transverse separation is 
of order Ib- The numerical coefficient 30 corresponds to 
Ib/ Pe ~ 10^^, a typical value for a galaxy cluster (assum- 
ing is a fraction of the radius). Being a logarithmic 
factor, the numerical value is insensitive to details. When 
electrons have moved a distance Lrr along the tangled 
field line, their three-dimensional root mean square (rms) 
displacement is i?*, where i?J ^ Lr^Ib ^ 301%. Beyond 
i?* , the motion of an electron is isotropic and uncorrelated 
with its previous path. 

Let us define as the time it takes for electrons to dif- 
fuse a distance Lrh along the field: ~ ^rr/i^W- ^'^^ 
t < t*, electrons diffuse anisotropically: I ~ (K||t)^/^,r 

Peexp[{K\\ty^^ /Ib]- For t > t^, however, electrons diffuse 
isotropically, and move in three dimensions according to 
R ~ where 



Rl/U - ilB/LRR)K\\ - lO-^KSp. 



(2) 



We see that for R > R^, the conduction is many orders 
of magnitude more efficient than when the field is ordered. 
However, Ki^ IS still a factor ^ 100 less than Ksp. The con- 
duction time is correspondingly 100 times longer than 
the Spitzer time tgp- As Table 1 shows, such weak con- 
duction is unlikely to have an important effect on cooling 
flows. 

The estimate given in equation is valid so long as 
X < Ib- This condition, which is likely to be satisfied by 
the gas in clusters (compare A with Ib ^ i?/(few) in Table 
1), ensures that collisions enable electrons to pass through 
magnetic mirrors caused by inhomogeneities in the field. 
If A > / B 1 conduction is supp r essed by an additional factor 
e <1 dChandran ct al. 1999| ; [Malyshkin fc Kulsrud 200l| ) 
since only a fraction of the electrons are able to penetrate 
the mirrors. This would cause the conduction time to in- 
crease by a factor 1/6. 



2.2. Conduction in a Multiscale Chaotic Magnetic Field 

A key assumption of the RR theory is the presence of a 
single Lyapunov length scale ^Lyap Ib- However, if the 
medium is turbulent, chaotic fluctuations will be present 
over a wide range of length scales. We generalize the the- 
ory for such a multiscale medium. 

We begin by re-expressing the single-scale theory as fol- 
lows. When two fleld lines are separated by a distance 
r smaller than Ib, their mean square separation (r^) in- 
creases with distance I along the field line according to a 
Lyapunov-like scaling. However, when r > Ib, the increase 
is given by the usual diffusion law, where (r^) increases by 
Ar^ ^ for a parallel displacement AZ '-^ Ib- Thus, we 
may describe the evolution in the two regimes, (r^) < 1% 
and (r^) > Z^, by the following two differential equations: 



d{r^)/dl - 2{r^)/h 



d{r^)/dl ~ II/Ib = h 



(3) 



Consider now a tangled magnetic field with a range 
of scales, and assume that the statistics of the magnetic 
field fluctuations are described by the Goldreich & Sridhar 
(1995, GS) theory of Alfvenic MHD turbulence. In a GS 
turbulent cascade, there is a range of scales l±_ perpendic- 
ular to the fleld, extending from a minimum scale Z_L,min to 
a maximum scale Ib- The fluctuations are anisotropic, so 
that for a given perpendicular scale Zj_ the corresponding 
parallel coherence scale l\\ is given by 



(4) 



For simplicity, we have selected the normalization such 
that on the outer scale Ib, the fluctuations are isotropic: 
l± In ^ Ib- The index a is equal to 2/3 for strong MHD 



turbulence (GS) and 3/ 4 for intermediate turbulence (Gol 
cich fc Sridhar 1997[ ). 

Using the single-scale equations (^ as a guide, it is 
straightforward to write down corresponding equations for 
a medium with a spectrum of fluctuations. We then iden- 
tify three regimes for the evolution of (r^). 

First, when r = (r^)^/^ < Z±,min, all the fluctuation 
scales in the medium contribute to Lyapunov-like growth 
(assuming they all behave chaotically), so 



d{r^) 
dl 



i/ii 



i/h 



din kj 
h 



2(r^ 

^ II .m in 



(5) 



where we have written k±_ = l/l±, and we have ignored 
constants of order unity in the normalization of the in- 
tegral. We see that the effective Lyapunov scale for the 
growth of r is the parallel coherence length of the smallest 
scale fluctuations in the medium, i.e. 



r ~ ro exp(Z//|i ,„i„). In 



'_L, mill's 



(6) 



Since usually Z||,niin Ib, the growth is rapid. 

Once r exceeds Z_L,min, the evolution switches to a second 
regime. We continue to have Lyapunov-like growth from 
scales l± > r, but there is diffusion-like growth for scales 
l±^ < r. The evolution equation for 
thus becomes 



d{r') ^ 2(^2^ dink, 
dl Ji/Ib h\ 



i/ii 

l/r 



/2 

dlnfci-t. (7) 
hi 
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If 3> r ^ ^±,min, each integral in equation is domi- 
nated by the scale r. Substituting equation in equation 
(0), we then obtain the solution 



r/lB ^ (I/Ib)'^". 



(8) 



Remarkably, the separation between two neighboring field 
lines becomes of order I b for a parallel translation of only 
Ib', this is much faster than in the RR theory which re- 
quires a parallel translation ~ 301b- The solution r oc 1^^° 
corresponds exactly to l± oc ^||^" in the turbulence model 
M). Thus, the rms separation of field lines grows along the 
'Aioldreich-Sridhar cone." 

When r > Ib, we enter a third regime, which corre- 
sponds to isotropic diffusion. 

^From equation (^, it is clear that thermal conduction 
in a multiscale chaotic field is almost as efficient as in 
Spitzer's theory. Replacing Lfm by Is in equation (^, 
we estimate Ktuib ~ (but see §2.3 for a better es- 



t imat I T' of the coe^fncient). As in ilic picvious subsection, 
we h ak^c aoQumod X < Is and have not included a mag 
netic mirror factor 6. In GS turbulence, perturbations 
on length scales l± < Ib have weak magnetic field fluc- 
t uat i c i us, AD/D {Ii/IbY^" < 1, ■lud cause 

neglig ible mirroring. — Only porturbationo on the aoalc Ib 



strong mirrori ng, but these have a n e gligible effect 



Kulsrid 2001) 



so lo ng go A < Ib -( Chandran ct al. lOOO ^- Malyohkin 



2.3. Numerical Solutions 

We have numerically integrated the differential equa- 
tions (|^) and (^, starting with an initial separation r = 
10^^'^Ib ~ Pe, and assuming a = 2/3 as appropriate for 
strong turbulence in the GS model. Figure 1 shows four 
numerical solutions for the evolution of r as a function of 
distance / along the field line, corresponding to four choices 
of the minimum scale Z_L.min of the turbulence. We see ex- 
ponential growth of r for I < Z^.min, and power-law growth 
for larger separations, confirming the scalings given in ^ 
and (|[). 

Let us define the decorrelation length Ldoc as the dis- 
tance along the field line for which the transverse sepa- 
ration r becomes equal to Ib- Figure 2 shows how Ldcc 
depends on Z^^min- When Z^^min ~ Ib, the turbulence is 
dominated by a single (outer) scale. This corresponds to 
the RR theory, and in this limit the decorrelation length is 



Thermal conductivity in a homogeneous magnetic field 
is known to be highly anisotropic — it is Spitzer along the 
field, but extraordinarily reduced in the transverse direc- 
tion. RR came up with the important insight that when 
the magnetic field is tangled and chaotic, thermal con- 
duction is enhanced significantly by the exponential di- 
vergence of neighboring field lines. However, even with 
this effect, CC estimated that the conductivity in galaxy 
clusters is below Spitzer by a factor > 100. 

We have shown in this Letter that if the field is chaotic 
over a wide range of length scales (factor of 100 or more), 
as might happen with MHD turbulence (GS), thermal con- 
duction is boosted to within a factor ~ 5 of the Spitzer 
value. Such strong conduction will have a significant effect 
on galaxy clusters (compare columns 6 and 8 in Table 1). 
It can transport heat to the center of a cluster to replace 
the energy lost through cooling, and it can also eliminate 
any thermal instability in the cooling gas. Thus, it may 
well reduc e the need for large-scale mass dropout in cool- 
ing fi ows (Tucker & Rosner 1983; Bertschinger & Meiksir 
1986). Only in the inner regions of some clusters (e.g 



3C295, Table 1) might there be significant dropout. It is 
worth noting that some authors have discussed potential 
problems wi t h invoking such strong co nduction ( p3inney"& 



Cowie 1981; Bregman fc David 1988), while others have 



2001) 



rebutted these arguments (Rosner ^ Tucker 1989; Bantoi 



An important requirement for the validity of our anal- 
ysis is that the magnetic field should behave chaotically, 
i.e., it should exhibit Lyapunov-like behavior over a wide 
range of scales. Weak MHD turbulence consists of a su- 
perposition of Alfven waves and is not chaotic. However, 
the model of strong and intermediate MHD turbulence de- 
veloped by GS is chaotic, as indicated by the breakdown 
of perturbation theory (Goldreich & Sridhar 1997). 

It should be noted that our theory of thermal conduc- 
tion does not require ongoing dynamic turbulence. Each 
episode of dynamic turbulence will leave behind a substan- 
tial level of frozen-in tangled fields even after the turbulent 
motions have ceased. Such leftover tangled fields should be 
sufficient to enhance conduction to the levels we estimate. 

As a final point, we should discuss a serious caveat. 
Chandra has found evidence for s harp temperature jumps 
in a few clusters, e.g., AbeU 2142 (|Markevitch et al. 2000D 
and AbeU 3667 ( |Vikhfinin et al. 2001| ). The observations 
indicate that con duction across the temperature jumps is 
far be low Spitzer (Markevitch et al. 2000; Ettori & Fabian 



that Ldcc is quite small, asymptoting to ^ 1.6/^. Since 
Ldec is the analog of the Rechester-Rosenbluth length Lj^ji 
for a multiscale medium, we may replace Lnn by Ldoc in 
equation (H) to estimate the diffusion constant in a turbu- 
lent medium'* 

Kturb ^ (iB/Ldcc) K|l ~ KSp/5, tturb ^ 5tsp. (9) 

Thus, if turbulence extends over at least two decades 
in scale, conduction is very efficient and approaches the 
Spitzer level to within a factor of a few. 

3. DISCUSSION 



large, las expected. However, for /^,min <; 10 ^Ib, we find 2000 ), in ap parent contradiction ^ with the estimate pre 



sented here. Vikhlinin et al. (2001) propose that the mag- 
netic field is stretched parallel to the interface, thus in- 
hibiting diffusion across the field on small scales. Extend- 
ing this idea, we suggest that when two distinct objects 
merge, as appears to be the case with the above clusters, 
the two regions (each of which is internally chaotic and 
highly conducting) may be thermally isolated from each 
other because their magnetic fields have not yet interpen- 
etrated each other. It is unclear how long such magnetic 
isolation will survive. 
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Fig. 1. — Variation of the rms transverse displacement r of two field lines vs longitudinal distance I along the field lines, for a = 2/3. 
Four choices of the minimum scale of turbulence are shown, ix,niin/'s = 10~^ (solid line), 10"^ (dashed), 10"" (~ pp//s, dotted), lO-^^ 

(~ Pe/^Bi long-dashed), which correspond to U^rainl^B = 2.2 X 10~^, 4.6 X 10"'', 4.6 X 10^*, 2.2 X 10~®, respectively. 



Fig. 2. — Variation of the decorrelation length Ldec vs Zj_ , 



Note that Ldec < 2Zb for , 



Table 1 

Comparison of the Cooling Time and the Conduction Time in Hydra A and 3C295 



Cluster R n kT X tcooi tsp Uurh 

name (kpc) (cm~^) (keV) (kpc) (Gyr) (Gyr) (Gyr) 



Hydra A 


100 


0.005 


3.6 


0.8 


5 


0.5 


2 


II 


10 


0.06 


3.1 


0.05 


0.5 


0.09 


0.4 


3C295 


100 


0.008 


5.0 


0.9 


7 


0.3 


2 


II 


10 


0.15 


3.0 


0.02 


0.3 


0.2 


1 



